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Fast evaluation and locality of overlap fermions 
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In order to construct improved overlap fermions, we start from a short ranged approximate Ginsparg- Wilson 
fermion and insert it into the overlap formula. We show that its polynomial evaluation is accelerated considerably 
compared to the standard Neuberger fermion. In addition the degree of locality is strongly improved. 



1. IMPROVED OVERLAP FERMIONS 

The Ginsparg- Wilson relation (GWR) for a lat- 
tice Dirac operator D reads 0| 

{D,75} = , [R local, {i?,75} ^ 0). (1) 

We choose R^^^y = j^^x.y, (a* > 0); ^-iid — 
75-D75, hence the GWR simplifies to 



A:=D- 11 



(2) 



It guarantees an exact — though lattice modified 
— chiral symmetry in even dimensions ||^, and 
parity symmetry in odd dimensions In both 
cases, the anomaly arises from the measure. 

Given some Dirac operator Z?o (local, no dou- 
blers), we can enforce the GWR on = -Dq ~ M 

by the overlap formula Aqv — /i^o / \J ^0^0 > so 
that Dov = Aov + A* obeys the GWR (for ^ in 
some allowed interval). The standard Neuberger 
fermion uses the Wilson operator Dq — Dyy @| • 

The guide-line for our concept to improve over- 
lap fermions is the following observation: if Dq 
is already a GW fermion (with respect to a fixed 
kernel R), then Dov — Dq- We can explicitly 
construct an approximate solution to the GWR, 
use it as Dq, and then we expect that the overlap 
formula only causes a small "chiral correction" , 



D, 



Do 



(3) 



Therefore, we first construct a short-ranged ap- 
proximate GW fermion, which includes elements 
of a truncated perfect action, and then we insert 
it as Do into the overlap formula. The result- 
ing Dov has an exact chiral symmetry, and we 

*Talk presented by W.B. at LATTICE 2001. 



expect a high level of locality, good scaling and 
approximate rotation invariance, because these 
properties hold for Dq, and we rely on relation 
(|^). Moreover, we predict a fast convergence un- 
der iterative evaluation of Dqu, because we start 
off in the right vicinity. 

All these properties have been tested and con- 
firmed in the 2-fiavour Schwinger model 
Here we report on convergence and locality in 
quenched QCD. 

2. THE HYPERCUBE FERMION (HF) 

Our construction of a short-ranged approxi- 
mate GW fermion is based on the perfect free 
fermion 0|, which is truncated to a unit hyper- 
cube. It consists of a vector term and a scalar 
term, D{x,y) = p^(a; - y)^^, + X{x - y), where 
Pi^{x ~ y)' Mx - y) ^ only for \xi - yi\ < 1, 
{i — 1 . . .4). The couplings are given in Ref. 
The first step is a "minimal gauging" of this HF, 
which means that the free couplings are attached 
to the shortest lattice paths only. This formu- 
lation provides already an excellent rotation be- 
haviour, but it suffers from strong additive mass 
renormalization 

To cure this problem, we attach a "link ampli- 
fication factor" 1/u to each link and tune it to 
Ucrit < 1 to approach at the chiral limit. Using 
different factors for the vector term and the scalar 
term further improves the GW approximation . 

Next we add a fat link, i.e. each link variable 
is substituted by a linear combination: (1 — a) 
"direct link" + a/6 "staple terms". Finally 
we also include a clover term. 

With this modest set of parameters we optimize 
the QCD spectrum at /3 = 6. In Ref. M we show 



2 



the full spectrum for typical configurations on a 
4"* lattice, as well as the eigenvalues with the least 
real parts on a 8"* lattice. 

3. POLYNOMIAL OVERLAP EVALUA- 
TION AND LOCALITY 

We consider two ways to evaluate the overlap 
Dirac operator by means of polynomial expan- 
sions, which have appeared in the literature. 
1) We introduce the Hermitean operator Hq = 
75^0, and evaluate 



Dov = A* 1 +7.-; 




(4) 



The last term corresponds to a sign function 
e{x) = 1 (—1) for X > (a; < 0), which is ap- 
proximated by a polynomial. 
2) We simply expand l/\/x by a polynomial, and 
insert AJAq for x. 

For the sign expansion, the convergence in the 
degree of the polynomial is exponential, if x G 
[— 1, 1]. Hence we first re-scale Hq Hq, so that 
the spectrum a{Ho) C [—1,1]. What matters is 
the density of eigenvalues (EVs) of Hq around 0, 
where the polynomial converges most painfully 
due to the discontinuity. For the HF, <t{Hq) has 
two peaks around ±1, which are shifted slightly 
towards the origin for Hq, leaving typically (at 
/3 = 6) a wide gap without EVs of Hq between 
about [—0.3,0.3]. On the other hand, Ho,w has a 
broad spectrum, ranging from about [—6,6], and 
after re-scaling one does obtains EVs near 0. His- 
tograms are shown in Ref. [^. 

For the l/\/x expansion the situation is simi- 
lar: here we re-scale so that x € [S, 1], {S > 0). It 
is crucial that 6 (the lower bound of the re-scaled 
spectrum) does not become very small; getting 
close to the singularity slows down the conver- 
gence. Again, the Wilson spectrum is very broad 
to start with, and re-scaling usually leads to small 
S, whereas it is kept larger for the HF. 

In Fig. 1^ we illustrate the convergence rate 
in QCD at /3 = 6 using Chebyshev polynomi- 
als. As an example, we discuss the e(a;) ex- 
pansion on the 4'' lattice. For Dq = Dw, the 
maximal (mean) deviation of the EVs of -D 
from the unit circle in C (with center 1) behaves 
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Figure 1. The maximal deviation of an energy 
eigenvalue of approximate overlap Dirac opera- 
tors (after re-scaling) from the unit circle in C. 



typically as d^""^ = exp(-0.134n) (d^-^"" = 
0.13 exp(— 0.134n)), where n is the degree of 
the Chebyshev polynomials. The correspond- 
ing decays for Dq = Dhf amount to c?^^"^ = 
exp(-0.737n) (d^l?" = 0.1 exp(-0.737n)). 

If we fix some affordable degree n, the precision 
of the GW approximation is superior for the HF 
by a factor = exp(0.6n), which takes a 

very considerable magnitude for realistic numbers 
like n = 20 . . . 100. 

On the other hand, we could fix a certain ac- 
curacy which we consider necessary, so the 
required degree n compares as nyy/nnp = 5-5. 
Referring to this number, the computational over- 
head of the HF, which amounts to a factor < 20, 
is compensated in part. It now takes only a rel- 
atively small progress in the scaling behaviour to 
make up for the remaining overhead. 

Fig. H shows that locality is clearly improved 
for the overlap-HF compared to the Neuberger 
fcrmion. Following Ref. [0, we measure the 
"maximal correlation" /(r) over a distance r. 
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Figure 2. Comparison of the degree of locality for 
different overlap fermions. 

We stay on the 12'' lattice and we observe also 
here a clear progress in the speed of convergence. 
This is related to the condition number c, i.e. the 
ratio of the upper and the lower bound of the 
spectrum of AJAq. Fig. ^ (on top) shows that the 
improved condition number for the HF is essen- 
tially due to the decrease of the upper bound. To 
illustrate the effect of c on the convergence rate, 
Fig. ^ (below) shows the deviation of /(r = 24) 
from the precise result if we use a moderate poly- 
nomial of degree n — 60. This deviation depends 
smoothly on c, up to an extra gap which makes 
the deviation for Dq = Dw yet a bit worse than c 
suggests, even at the optimal parameter /i — 1.64. 

4. CONCLUDING REMARKS 

Compared to the standard Neuberger fermion, 
the overlap-HF gains significantly in the conver- 
gence rate under polynomial evaluation. We have 
studied the evolution of the full operator in the 
polynomial at /3 = 6, which might be roughly 
comparable to the situation at stronger coupling, 
say P — 5.75, after treating the 0(10) worst 
modes separately, as it is often done. 

Moreover, we observed a superior degree of lo- 
cality. This suggests that the overlap formula is 
applicable up to stronger coupling (the limit for 
P was discussed in Ref. Q). 

The scaling quality remains to be tested; we are 
currently measuring meson dispersions. Another 
important question is the applicability of precon- 
ditioning techniques. They were already applied 
successfully to the HF without overlap ]11|. 
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Figure 3. On top: history of the spectral bounds. 
Their ratio ( — condition number) is essential for 
the polynomial convergence rate (see below). 
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